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Spectral energy distributions and high-energy emission of 
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ABSTRACT 

Based on identifications from the Veron and Quasars.org catalogs, we determine the 
optical-to- X-ray spectral indices for a sample of 201 BL Lac type objects (BLLs) and 
search for trends in the distribution of these indices of the sources detected in high- 
energy bands. We find that EGRET-detected sources are low-energy peaked and that 
the positional correlation with the arrival directions of ultra-high-energy cosmic rays 
from the previously studied AGASA, Yakutsk and High Resolution Fly's Eye samples 
is dominated by low-energy-peaked BLLs. 



Key words: galaxies: BL Lacertae objets: general 
tion of particles. 
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1 INTRODUCTION 

BL Lac type objects (hereafter BLLs) attract consider- 
able attention of both astrophysicists and particle physi- 
cists. Among other reasons, it is caused by observation 
of vcry-high-en ergy (~ 1 TeV) g amma rays fr om these 



sources (see e.g. lAharonia ^ j200dhlWagnei1 (20071 ) and ref- 
erences therein) and by claims for potential association of 
BLLs w ith ultra-high-energy ( UHE ; > 10 1 9 eV) cosmic rays 
(CRs) fiTnvakov fc Tkachevl l200ll . 120021 : iGorbunov et al.l 
l2002] , l2004l ; lAbbasi et alj2006h . Most probably, only particu- 
lar subclasses of BLLs are responsible for these high-energy 
emissions. More precise determination of these subclasses 
may shed light on the intrinsic mechanisms of particle accel- 
eration in blazars and may help in determination of potential 
sources yet unobserved at high or ultra-high energies. 

Broadband spectral energy distributions (SEDs) of 
BLLs are known to have a very specific non-thermal 
two-bump shape. In a popular model, the two bumps 
are caused by the synchrotron and inverse-Com pt on ra - 



diations (see e.g. Maraschi. Ghisellini fc Celottil (|l992l) 



ISikora. Begelman fe Reed i| 19941) ). Their origin i s attributed 



to oth er mechanisms in some models (see e.g. iMannheiml 
|l993t )). While the overall two-bump shape is quite com- 
mon, the location of the peaks varies strongly from source 
to source. In particular, the well-measured low-energy peak 
corresponds to infrared or optical frequencies in the so-called 
low-energy-peaked BLLs (LBLs) and to X-ray frequencies in 
the high-energy-peaked objects (HBLs). As it has been re- 
cently understood, intermediate cases are also present. The 
correlation of the peak positio n with the intrinsic power of 
the object l|Fossati et al.l ll998) is currently under discussion 
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(see e.g. iPadovanH (|2007l )). The second, often worse mea- 
sured peak is located in the gamma-ray band (MeV to GeV 
in LBLs and hundreds of GeV in HBLs). 



Since LBLs are bright in the optical band and relatively 
faint in X rays (keV frequencies correspond to a dip between 
two peaks in this case) while HBLs peak in the X-ray band, 
the optical-to- X-ray broadband spectral index aox provides 
a good measure of the posit ion of the first peak in SED 
(see e.g. iDonato et al.l ([2001 ) ) . Detailed measurements of 
the SED are availabie for a small fraction of BLLs while to 
know aox one needs only optical and X-ray measurements 
performed for a much wider sample of sources. In this note, 
we use cvox as the quantity which characterizes the SED of 
a BLL; results of a more detailed multiwavelength study will 
be reported elsewhere. We determine aox for a large sample 
of BLLs and look for a correlation between the index and the 
high-energy emissivity of a source. We will see that EGRET- 
detected sources are mostly LBLs. One of the primary goals 
of this study is to specify the class of BLLs which correlate 
with the arrival directions of ultra-high-energy cosmic rays; 
we will see that the correlations observed previously are also 
saturated by LBLs. 



The rest of the Letter is organized as follows. In Sec. [2j 
the sample of BLLs is discussed and the index aox is de- 
termined. In Sec. [3] we demonstrate the trends in the dis- 
tribution of aox of EGRET sources. In Sec. [4] we briefly 
review previously reported UHECR - BL Lac correlation 
and discuss the distribution of aox of correlated objects. 
Brief conclusions are summarized in Sec. [5] 
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2 THE SAMPLE 

We stu dy objects classified as c onfirm ed BLLs (class BL or 
HP) in IVeron-Cettv fc Veron! |2003l ) (hereafter the Veron 
catalog). Objects from this catalog were searched for X- 
ray i dentifications in ROSAT d ata in the Quasars.org cat- 
alog (|Flesch fc Hardcastldl2004h . There, the Veron objects 
have been cross-correla t ed wi t h vari ou s publish e d ROS AT 
catalogs JVoges et all ll 19991. l2000l l: iROSAll (|2000bl lah: 
IWhite. Giommi fc Angelinil (|2000l )f and the probability of 
the true identification was presented for each case. Of the 
sources from different catalogs associated with a given BLL, 
we select one with the highest probability; we drop the 
source from the sample if this probability is below 68%. In 
such a way, we obtain a sample of 201 objects with V-band 
magnitudes given in the Veron catalog and X-ray identifi- 
cations with one of the ROSAT catalogs. We note that in 
some cases, the name of the object in the Veron catalog has 
the ROSAT prefix but this identification has very low prob- 
ability according to Quasars.org and the corresponding BLL 
was therefore dropped from the sample. 

For t he objects identified with the ROSAT la test PSPC 
catalogs jVoges et al.ll 19991 . 120001 ; [ROSATll2000blL we calcu- 
late the flux in (0.1 -r 2.4) keV band from the coun t rate 
and the hardness ratio following IVoges et al.l (|l999h . In a 
few cases when the best ide n tification is with the catalog of 
IWhite. Giommi fc Angelinil l|l994 l2000t ). we take the flux 
value presented there. If the object was observed by the 
high- resolution imager ( HRI), th e best i dentification is of- 
ten with the catalog of iROSATl (|2000al); we use the flux 
values given in the BMW catalog ( Panzera et al.l 120031 ) in 
these cases. These latter fluxes correspond to the same en- 
ergy band and assume a Crab-like spectrum for each source. 
It is a relatively rough approximation caused by a poor spec- 
tral sensitivity of HRI; however, for firm identifications, the 
BMW and PSPC fluxes are well correlated. 

For each object in the sample, we calculate the optical- 
to-X-ray spectral index «ox using the textbook definition 
l|Carroll fc Ostliell2007l ), 



Fo 
Fx 



where Ft dvi is the amount of energy with frequencies be- 
tween Vi and Vi + dvi per unit area per second, observed 
by a detector aimed at the source and measured e.g. in 
W/(cm 2 -s); index i stands here for either l/-band, i =0, 
or X-ray band, i =X, frequencies and fluxes. 

Most of the BLLs are strongly variable at all wave- 
lengths and therefore the use of non- simultaneous observa- 
tions may introduce random errors as large as 0.3 to aox- 
Given the fact that the X-ray sample is essentially flux- 
limited, these errors may be asymmetric for faint objects. 

An important correction to V and therefore to aox may 
arise from the contribution of the host galaxy emission. We 
use the magnitudes from the Veron catalog which are not 
corrected for the contribution of the host galaxies. These 
corrections might cause significant systematic errors in aox 
given the fact that the host galaxy contribution is more im- 
portant in the visible light than in X rays. The effect of this 
correction is stronger in fainter, less beamed or misaligned 
BLLs. Account of this correction on the case-by-case basis 
would result in a significant reduction of the statistics since 



for many objects in the sample (notably for those associated 
with EGRET sources or correlated with UHECR) the host 
galaxy was not detected, indicating either powerful or well 
beamed source. Another approach is to use infrared colour 
indices to estimate, on a statistical basis, whether th e con- 
tribut ion of a host galaxy is important. Developed by I Glass! 
|l98lf ). t his approach was system atically applied to Veron 
BLLs bv lChen, Fu fc Gaol (|2005l ). where J, H and K mag- 
nitudes were obtained from the 2MASS survey and the ob- 
jects for which a significant contribution of the host galaxy 
is expected were determined. To estimate the effect of the 
contamination by the emission of the host galaxies on our 
results, we removed these objects (37 out of 201) from our 
sample and checked that the effects we advocate are present 
in the reduced sample as welQ, see below. Since the starlight 
of the host galax y is expected to peak in the i nfrare d (see 
e.g. discussion in iKotilainen. Falomo fc Scarpa! (|l998h ). the 
effect in the IZ-band should be even smaller. 

Our sample is derived from the Veron catalog which 
is incomplete. Therefore, the results of the study may be 
biased; it is not guaranteed that they are generic and hold for 
all BLLs in the Universe. However, by comparing different 
subsets of one and the same catalog, we partially remove 
the effects of the incompletness and reveal important trends 
which can be tested in future studies. Though the selection 
effects may be different for different subsamples, the most 
important of them are under control in the present study. 



3 GAMMA-RAY EMITTERS 

The identification of BLL counterparts of the EGRET 
sources is a nontrivial task and becomes sometimes a 
subject of debates. There exist EGRET sources iden- 
tified with BLLs by one authors but identified with 
other objects (e.g. clusters of galaxies) by others. For 
this study, we selec t all p otential BLL identifications 
from Hartman ct al. (1999); iMattox. Hartman fc Reimerl 
|200ll ); [Go"rbunov et al.l (|2002l ) and divided them into "high- 
probability" and "others" following the id e ntific ation prob- 
abilities from lMattox. Hartman fc Reimerl (|200ll ) (the high- 
pr obability subsample includes object s from Tables 1 and 2 
of lMattox. Hartman fc Reimerl (|200l| j). 

In Fig. [T] we present magnitudes and spectral indices 
of those BLLs in our sample which are associated with 
catalogued EGRET sources. Fig. [2] compares the distri- 
bution of the spectral indices of correlated objects with 
that of all objects in the sample. In agreement with the 
two-bump synchrotron-self-Compton SED model, EGRET- 
selected sources are low-energy peaked. To quantify this 
statement, we compared the distribution of aox for the 
EGRET BLLs and that for all objects in the sample by 
means of the Kolmogorov-Smirnov (KS) test. The probabil- 
ity that the parent distribution is the same is P ~ 5 ■ 10~ 3 ; 



it reduces to ~ 10 if only high-probability identifications 
are used or if the BLLs with potentially strong contribution 



1 Of our sample , 34 o bjects are not included in the catalog of 
IChen. Fu fc Gaol d2005l> . We assume that the lack of infrared iden- 
tification implies negligible contribution of the host galaxy for 
these sources. 
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Figure 1. Optical-to-X-ray spectral index aox versus V-band 
magnitude for all BLLs in the sample (gray) and for EGRET- 
detected objects (black). Black boxes correspond to prob- 
able identifications (according to Mat tox. Hartman fe Reimerl 
teOOlfO, triangles correspond to le ss probable or not studied by 
iMattox. Hartman fc Reimerl J200lf) objects. 
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Figure 3. Optical-to-X-ray spectral index aox versus V-band 
magnitude for all BLLs in the sample (gray) and for those cor- 
related with arrival directions of ultra-high-energy cosmic rays 
(black). Black boxes denote objects correlated with AGASA and 
Yakutsk cosmic rays (large boxes correspond to objects correlated 
with doublets), triangles denote objects correlated with HiRes 
cosmic rays (see text for more details). 
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Figure 2. Distribution of «ox f° r al l BLLs in the sample (gray) 
and for objects possibly detected by EGRET (dashed). 



of the host galaxy are dropped. The change in the KS prob- 
ability corresponds to the decrease of the sample size. We 
note that all 14 high-probability EGRET identifications cor- 
respond to BLLs without significant expected contribution 
from the host galaxy, so a more detailed account of this bias, 
would it be possible, may only strengthen our conclusion. 

A similar study for the BLLs detected at very high en- 
ergies (E > 200 GeV) is hardly possible because the^o are 
not drawn from a full-sky survey but from pointed observa- 
tions of X-ray selected objects (low aox). One should note 
however that strong X-ray flux does not guarantee that a 
BLL is a TeV emitter; indeed MAGIC performed a dedi- 
cated search for TeV emission from X- ray selected blazar s 
and in many cases did not discover it (I Albert et al.l [20071 ). 
We note in passing that no significant preference in aox is 
found for the BLLs which are possib le emitters of 10-GeV 
gamma-rays l|Gorbunov et al.l l2005al ); they have, on aver- 
age, intermediate values of aox between EGRET and TeV 
emitters. 



2 The list of these sources may be found on the MAGIC webpage 
http : //www. mppmu. mpg. de/~ rwagner/sources/index. html 



4 COSMIC-RAY CORRELATION 

Previous studies reported significant correlation between 
various samples of BLLs from the Veron catalog and 
various samples of UHECRs. Some of the studies 
l|Tinvakov fc Tkachevl [20021 ; IGorbunov et alj|2002f ) used re- 
construction of the arrival directions in the Galactic mag- 
netic field assuming charged cosmic-ray particles; we will not 
discuss them here be cause of th e ambi guity in the magnetic- 
field models (see e.g. ICao et ail l|2007h ). On the other hand, 
a number of studies suggest correlation which assumes zero 
deflection (neutral particles); these results are less model- 
dependent and much more intriguing because neutral UHE 
particles from BLLs would challenge conventional models 
of cosmic-ray physics. These claims include the correlation 
found in the sample of cosmic rays observed by the Akeno 
Giant Air- Shower Array (AGASA; sample with estimated 
primary energies E > 4.8 • 10 19 eV) and the Yakutsk Ex- 
tensive Air Shower Array (Yakutsk; E > 2.4 ■ 10 19 eV) de- 
tectors where an excess of pairs 'BL L - cosmic ray' was 
seen at separations less than 2.5° (jTinvakov fc Tkachevl 
l2001f l and in a sample of events with E > 10 19 eV ob- 
served by the High Resolution Fly's Eye detector (Hi Res) 
for separations less than 0.8° l)Gorbunov et al.l 120041 ). In 
both cases the separation was consistent with the detec- 
tor's angular resolution (which was much better in HiRes 
than in AGASA and Yakutsk). The correlation with the 
HiRes sample was confirmed in an unbinned study and 
was f ound to be held at lower energies using unpublished 
data ()Abbasi et al.l 120061 ). The probability to observe the 
correlation with t hree independent expe r iment s by chance 
was estimated by IGorbunov fc Troitskvl (|2005r ) as 3 ■ 10 -5 
by a Monte-Carlo study which took into account statisti- 
cal penalty for multiple tries (various catalogs of potential 
sources tested for correlation). 

In Fig. [3] we present magnitudes and spectral indices of 
those BLLs in our sample which are located within 2.5° of 
the arrival directio ns of AGASA and Yak u tsk co smic rays of 
the sample used in iTinvakov fc Tkachevl (|200ll) and within 
0.8° o f those of the HiRes sample used in Gorbunov et al.l 
l|2004l ). Fig. [4] compares the distribution of the spectral in- 
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Figure 4. Distribution of aox f° r all BLLs in the sample 
(gray) and for objects correlated with ultra-high-energy cosmic 
rays (dashed). 



dices of correlated objects with that of all objects in the 
sample. We see that the correlation is dominated by low- 
energy peaked BLLs; the KS test for the two distributions 
gives P ftf 4 • 10 -3 . This co nfirms and improves ou r pre- 
liminary approximate results l|Gorbunov et al.ll2005bl ) based 
on a sample of BLLs with X-ray identifications 'by eye'. 
The results of the present study (cf. Fig. [2] and Fig. [4} 
support also the asso ciation between UHEC R and EGRET 
sources suggested by iGorbunov etaH l|2002t ). This associa- 
tion is not surprising since both acceleration and propaga- 
tion of UHECR are inevitably accompanied by emission of 
secondary energetic photons which in turn interact with the 
cosmic background radiation and lose their energy until it 
reaches ~ 0.1 -r- 10 GeV. 

The discussion of the potential bias due to host galaxies 
is very similar to the case of the EGRET sources. Restriction 
of the sample to objects with no expected host galaxy con- 
tribution results in the KS probability of P w 10 -2 , change 
being consistent with the decrease of the sample size. Of 21 
objects in the sample associated with UHE CRs, only 3 may 
be affe cted by the host galaxy, according to lChen. Fu fc Gad 
20051). This fact is in agree ment with previous observations 



Tinvakov fc Tkachevll200"lh that the UHECR correlation 



stronger for objects with unknown redshifts. 

The co rrelation between BLL s and UHECRs seen in 
HiRes data (jGorbunov et al.l|20 04f) has been tested recentl y 
by the Pierre Auger (PA) Collaboration (jHarari et al.ll2007h : 
no positive signal was found. This is not conclusive how- 
ever for the following reasons. Firstly, PA is located in the 
Southern hemisphere and sees different BLLs than other 
experiments; moreover, due to incompleteness of the cat- 
alogs, the number of potential UHECR emitters is much 
less in the South. For instance, of 99 objects in our sample 
which have qox > 1 (which seem to correlate with UHECR 
stronger, see Fig. [4}, 82 can be seen by HiRes but only 46 
are in the field of view of PA, most of them only at large 
zenith angles and for a small fraction of time. The 'factor 
of merit' depends also on the angular resolution of the ex- 
periment which is twice worse in the PA surface detector 
than in HiRes (stereoscop i c mod e). This problem was quan- 
tified by IGorbunov et al.l (|2006t ) where it has been shown 
that to reach th e HiR es sensitivity for the signal found by 
IGorbunov etafl (|2004l ) in the set of 271 events, PA has to 
accumulate ~ 3500 events in the same energy range. We 



note that lHarari et al. (2007) used 1672 events for the test 

of this signal. 

Second ly, as it has be e n poin ted out bv lGorbunov et all 

d2004h and febbasi et all d2006h and further discussed by 
iTinvakov fc Tkachevl ^2006). the correlation observed by 
HiRes implies neutral cosmic particles travelling for cos- 
mological distances, the fact which requires unconventional 
physics. Most probably the primary particles of the re- 
sulting air showers are neither protons, nor nuclei. How- 
ever, the energy determination of the PA surface detec- 
tor is extremely sensitive to the type of the primary cos- 
mic particle because of very strong sensitivity of water 
tanks to muons in the air shower. For instance, energies 
of gamma ra ys are always underestimated by a factor of 
a few (see e.g. iBilloir. Roucelle fc Hamilton! |2007l )'). Due to 
the steeply falling spectrum of UHECRs, this may dilute 
the observed signal. PA possesses fluorescent detectors sim- 
ilar to those of HiRes, but the BLL correlation was never 
tested with them. Future tests of the correlation should be 
performed with high-resolution fluorescent detectors, prefer- 
abl y in the Northern hem isphere, the Telescope Array (see 
e.g. iMartens et al.l l|2007h ) providing a good example. 



5 CONCLUSIONS 

We have studied the distibution of broadband optical-to- 
X-ray spectral indices aox of 201 confirmed BLLs from the 
Veron catalog which have at least 68% confident X-ray iden- 
tification in the Quasars.org catalog. In accordance with the 
synchrotron-self-Compton two-bump SED models, a sub- 
sample of EGRET-detected BLLs has on average high aox 
(sources are strong in optical but relatively faint in X rays). 
One of the most important results of the study relates to the 
BLLs correlated with ultra-high-energy cosmic rays: most 
of them are low-energy peaked. Both for the EGRET and 
UHECR associations, the contribution of the host galaxy is 
expect ed to be negligible accor ding to the infrared photom- 
etry of IChen. Fu fc Gaol l|2005h . This fact agrees well with 
the expectations that the EGRET sample is dominated by 
beamed sources and that the UHECR sources are powerful 
and possess strongly collimated jets pointing precisely to the 
observer. The present study narrows the class of potential 
UHECR emiters to those with high aox. This fact may be 
used for further tests of BLL - cosmic ray correlation with 
fluorescent detectors and may shed light on the origin of 
the correlated cosmic-ray particles. If confirmed, this corre- 
lation would point to completely new phenomena in particle 
physics or/and astrophysics because no known neutral par- 
ticles of these energies are expected to travel for distances 
larger than ~ 10 Mpc. 
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